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Abstract 



K.^ ' We study charging effects and tunneling in the single electron box. Tun- 

Q"^ ' neling mixes different charge states and in the nonperturbative regime the 

o' 

t^^ ' charge in the island may be strongly screened. When charge states are nearly 



degenerate the screening of the charge is strong even in the weak tunneling 



C^ ' regime. Virtual tunneling processes reduce both the level splitting A and the 

a" 

tunneling strength a. The charge on the island and the decay rates are calcu- 



lated. In the strong tunneling regime also nondegenerate states are affected by 
tunneling. Strong-coupling scaling renormalizes the effective capacitance, a 



^ ' result which we confirm by Monte Carlo simulations. The tunneling strength 



a scales to smaller values into the regime where the weak-coupling scaling ap- 
plies. We propose a two stage scaling procedure providing the unified picture 
for the problem. The scaling analysis is also extended to superconducting 

tunnel junctions with finite subgap conductance. 
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I. INTRODUCTION 

Charging effects strongly modify tlie transport properties of systems of small tunnel junc- 
tions. The transfer of an electron requires a typical electrostatic energy cost Ec = e^/2C, 
so at sufficiently low temperatures (T ^ Ec) charge transfer is suppressed, a phenomenon 
known as the Coulomb blockade |[T|-§[. Tunneling induces quantum fluctuations of the charge 
in the electrodes. Strong fluctuations screen the charge even at T = and charging effects 
are weakened. Quantum fluctuations are strong if the typical tunneling resistance R is small 
enough (oq = 2n%R ^ ■'-)• Even in the weak tunneling regime Oq < 1; quantum fluctuations 
are strong if the lowest lying charge states are nearly degenerate. 

A simple and widely studied device where these effects are manifest is the the single 
electron box [|], shown in Fig. 1. It consists of a normal metal island connected via a 
capacitor and a tunnel junction to a voltage source. The energy gap A between the two 
lowest lying states can range between Ec and zero and is controlled by the external voltage 
(see Fig. 2). If the tunneling is strong {a^ > 1) the screening of the charge results in a 
reduction of the overall effective bandwidth Ec- Also for weaker tunneling, ao < 1, infrared 
divergent tunneling processes strongly mix nearly degenerate charge states. As a result the 
gap A is reduced near the degeneracy points, but the overall bandwidth Ec is unaffected. 
The quantitative analysis of these regimes lies beyond the scope of perturbation theory. The 
problem has recently received much attention P-§, however several previously published 
results i5|,|ToHl2| are in mutual disagreement. Our motivation is to provide a plausible unified 
description for this problem. 

In section II we introduce the model. In section III we discuss the scaling in the weak 
tunneling regime. In section IV we review the scaling in the strong tunneling regime and 
propose a two stage scaling procedure to provide a unifying picture for large and small aQ. 
We also obtained the strong coupling scaling of the effective capacitance by a Monte Carlo 
simulation. Then we compare with other non-perturbative techniques in the ao > 1 regime. 
In section V we present results for other observable quantities, we extend our analysis to 



finite temperatures and finally discuss the scaling for a superconducting junction with subgap 
quasiparticle tunneling. 

II. THE MODEL 

In the absence of tunneling the thermodynamics of the single electron box is governed 
by the electrostatic energy Eq{Q) = Q"^ /2C where C = Cj + Cs (see Fig. 1) and the charge 
Q = {Qx — ne) is composed by the n excess electrons in the island and the continuous 
"external charge" Q^ = CgV^ induced by the voltage source. The measurable voltage at 
the junction is {V) = {Q)/C. The energy spectrum of the system as a function of Q^ is 
shown in Fig. 2. The lowest lying levels are degenerate at Qx = {k + l/2)e. At T = a 
well defined n is selected, which depends on Q^ and changes by ±1 when Q^ crosses the 
degeneracy points. The result is a step structure of (n) and a sawtooth shape of {V) as a 
function of Qx 0- At finite temperature both are smoothened by fluctuations. 

In the presence of tunneling the single electron box can be described by the following 
Hamiltonian 

H = ^ + I] e(p) C+^Cp,, + J2^(^)4aCka, + Y.^kp4^Cp„ + h.C. . (1) 

po- ka kpa 

Here n = J^pa (^pa^po- 5 c^^ ^^^ (^ta ^^^ creation operators of the electrons in the island and 
in the lead respectively, and e{p) and e{k) are their kinetic energies. We consider a wide 
tunnel junction where tunneling takes place through A^ independent channels. In this case 
the tunneling strength is given by oq = N\T\'^p'^ , where T is the tunneling amplitude for a 
single channel and p is the densities of states per channel, which is assumed to be equal in 
both electrodes. 

In order to concentrate on the quantum dynamics of the charge we integrate out the elec- 
tronic degrees of freedom |^ in the partition function Z(Q§). This is achieved by decoupling 
the quartic charging term via a Hubbard- Stratonovich transformation which introduces the 
phase (p{t). After the integration ^(Q§) can be expressed as a path integral depending on 



the phase only. Expanding the effective action in powers of the tunnehng amphtudes Tkp we 
obtain a 1/A^ expansion, i.e. the 2n*^ order term is proportional to A^lTp'^p^" = Oq / N"^'^ . 
Since in our case the number of channels A^ is large and the nominal conductance ao is finite 
it is sufficient to retain only the first term of the expansion, which is given by the well known 
path-integral representation [|T^ 

2iQ.) = E e^'^^"^^^/^ / dcfo / V^ir) e-^t^^^" , (2) 

m=~OD ''-°° ''fo 



with the effective action 

S[^{t)] = l'dr^(^^^ - f^dr j' dr' a{r - r') cos [^(r) - (^(r')] • (3) 

The Fourier transform of the dissipative kernel is a{ijJn) = — 7rQ;o|ti;„| up to a high energy 
cutoff C(Jc. The summation over winding numbers, ^{{3) = f{0) + 2TTm , reflects the discrete- 
ness of the charge. The external charge Q^ can be viewed as a gauge field and appears in 
the phase factor together with the winding number m . 

It is useful to consider also the dual representation expressed in terms of charges |0] 

oo ^(j —ne '^^ 1 r3 3 

Here Q{t) = Qx—ne~ e Z)?, ©(t— Tj) — 0(t— r/) . Equation (4) maps the tunneling problem 
on a gas of interacting blips and antiblips (see Fig. 3) each one representing a tunneling 
transition at time Tj in which the charge changes in units of e. The tunneling events force 
a rearrangement of the other electrons in the two electrodes, which is represented by the 
a(r) lines. 

We stress that both equations (2-3) and (4) are exact representations of the original 
problem in the limit of interest, N ^ oo and ao finite. The extension to finite values of 
A^ has also been considered previously pO|. 



III. WEAK TUNNELING REGIME 

In the weak tunneling regime ao < 1 it is convenient to start from the charge rep- 
resentation Eq. (4). Due to the symmetries of the system we can focus on the interval 



< Qx < e/2 . The main controlling parameters are the strength of the tunneling ao and 
the energy difference between the two lowest charge states 

Ao{Qx)=EoiQx-e)-EoiQx) = Ecil-2Qx/e) . (5) 

Due to the infrared divergent behaviour of q;(co'), perturbation theory in oq breaks down near 
the degeneracy points. Indeed the leading contribution to the ground state energy, which 
up to singular terms is Eq = Eq{Q,j.) — oqAq ln(Ao/ci;c) shows an infrared singularity for 
Ao ^ 0. Hence both (V) = -^(c/EgV^^o) and (n) = (g^-C(1/))/e are logarithmically 
divergent. 

Close to the degeneracy points, Ao(Qx) ^ Ec the low energy physics is determined by 
the lowest two states and involves tunneling events with energy differences smaller than Eq- 
Thus the cutoff of the model, Uc, is approximately given by Eq- In this two-level-system 
(TLS) approximation only the terms n = 0,1 of Eq. (4) are retained Z = J2n=o,i Gn{P) and 
the only trajectories Q{t) contributing to Qn{f3) are such that blips and antiblips alternate. 
For instance we retain the diagrams (a), (6) and (c) in Fig. 3 and drop (d). The bare 
propagators are ^„(r) = exp(— |r| Eq{Qx — ne)). 

The perturbative corrections are logarithmic in the regime Tc ^ t. The first order term 
(blip-antiblip pair) in the leading logarithmic approximation reads as 



Z(i) (g., r) ^ e-^°(«^)^ (l - ao(Aor + 1) \n{u:,f) 

+ e-(^«(«-)+^«)^ (l - «o(-Aor + 1) ln(cu,f )) (6) 



where 1/f = max[l/r, Aq, T] acts as a low frequency cutoff which regularizes the infrared 
singularities and 0{1) constants depending on the details of the high frequency cutoff 
procedure have been ignored. 



Our model Eq. (4) is similar to the model used by Anderson et al. [|T^ to study the 
single-channel S = 1/2 Kondo problem. We could use the same renormalization group 
(RG) technique to treat the infrared singularities, namely progressively eliminating close 
blip-antiblip pairs by increasing of the short time cutoff Tc = l/wc H]. However, since 



there are technical differences (e.g. here the interaction is pairwise, and the short time 
cutoff procedure is chosen such as to guarantee that a{ujn = 0) = 0), some care is needed 
when using the scahng technique of Ref. |T^. We chose to ehminate high frequencies; so 
we split the kernel a(r) into slow and fast parts, the latter containing the frequencies we 
want to integrate out, a(r) = as(r) + ^/(t) . In first order in the tunneling strength we 
look at all the configurations containing only one "fast" a-line which can connect either 
a close blip-antiblip pair or a pair separated by some blip-antiblip insertion. Only the 
former are effective in the first step of the renormalization. Indeed the integration of the 
complete (fast and slow modes) line of the latter kind does not give infrared singularities 
because of phase space restrictions, thus we discard them. Interaction lines overarching 
blip-antiblip insertions in higher order diagrams (eg. the "rainbow" of Fig. 3b) will enter 
the renormalization in the successive steps, once the blip-antiblip insertions are eliminated. 
In this procedure the running parameters will never be renormalized by diagrams containing 
crossing a-lines. Our treatment is equivalent to the "Non-Crossing Approximation" and can 
be justified by direct calculation. Indeed elimination a close pair with crossing a-lines (see 
Fig. 3c) yields an interaction proportional to r^^ which does not lead to any logarithmic 
singularity. The bare ground state propagator Qo{t) is renormalized in lowest order by 
the fast modes Uc — Suc < uj < uo^ of a single blip-antiblip close pair which generate a 
contribution — ao(l + ^ot){^^c/^c) ■ A similar expression gives the renormalization of the 
propagator for the first excited state Qiij). 

The general n-th order terms can be cast in a way such that the partition function 
preserves the original form if the gap and the interaction scale as 

^^P^ = (2aZ^ - 1)^ . !^iM!) = (2«Z'f . (7) 

dmojc ^ ' ijJc dm ijJc ^ ' 

The renormalized ground state energy is given by Eq = Eq{Qx) + l/2(Ao — A) . Here Z 
is the wave function renormalization which enters always together with the interaction in 
the combination 2aZ'^ , and the —1 in the first equation enters because of the dimensional 
factor Ur . 



The scaling equations can be readily integrated down to a low energy scale Uc 

a{uj,) = ao (l + 2ao In (—))"', A(a;,) = Aq (l + 2ao In (—))"' . (8) 

At zero temperature the renormalized gap provides the low energy cutoff so we have to 
stop the scaling at lo^ = A. Then Eq. (8) becomes a self-consistent equation for A. At 
finite temperatures, if T > A the RG has to be stopped at Uc = 'max[A,T] , where the 
infrared singularities in the perturbation expansion disappear. 

Typical solutions are shown in Fig. 4. Notice that the renormalized gap A does not 
vanish. 

IV. STRONG COUPLING REGIME 

In the regime ao ^ 1 the charge fluctuates strongly so it may be convenient to start from 
the phase representation Eqs. (2-3). We analyze the problem by various nonperturbative 
techniques. We start with the RG analysis of Refs. P, p!5| , p!6| . This treatment is restricted to 
the 171 = sector of the partition function Eq. (2), so the discreteness of the charge is not 
explicitly accounted for. However, we will show that the overall conclusions of the scaling 
theory can be carried over to the general case, described by Eq. (2). 

The scaling equations are obtained perturbatively in 1/a p|, p^p!6| . 

^^i^sM = (i _ 1) (£,/,j ; §ilM._i/a^ , (9) 

a In ujc ^a '' a m Uc 

where a = 27r^a. In the present approach we determine the scaling of the overall bandwidth 
Eq, but not that of the gap A, exactly the opposite of the weak coupling region. However we 
assume that the gap A scales similarly to the bandwidth. Equation (9) shows that a decays 
in the strong coupling regime towards weak coupling; so in the case of strong tunneling a 
two step scaling procedure is called for which provides the desired unified picture. As the 
energy scale Uc (e.g. temperature) decreases, first a decays from its large initial value Oq 
according to Eq.(9) 



-^— = 11^ (10) 

a (uc) 1 - (1/ao) In (t^co/^c)) 

to a ~ 1. The bandwidth also decreases upon renormahzation and for a ~ 1 reaches the 
value E^ ^ E^ exp (—do) • At this point the cutoff reaches Uc ~ E^. From there on we 
use the weak coupling scaling of section III with initial values u^ ~ E^ , the exponentially 
suppressed bandwidth, Aq ~ Aq exp(— do) , and do ~ 1 implying oq ~ l/27r^. The flow 
of a (ujc) will be governed by Eq.(7). Again the tunneling strength decreases and eventually 
we stop the RG when we reach the low-energy cutoff. The final formula for the gap in the 
tto ^ 1 regime is 

Ao e""° 
^^'"^^= l + .r-Mn(EOe--oK) " ^''^ 

Next we checked the exponential suppression of the bandwidth by a Monte Carlo simu- 
lation starting from the phase representation in the normalized form 

1 °° r 

ZiQ,) = ——-- E ^Me^-"^^^/^ ; Z{m) = V6^iT)e-^-^'^^ (12) 

where we decomposed ^{j) = Sipi^r) + ipo + 27rmr//? , Sm [^fiT)] = <S[Sf{T) + 27rmr//?] . 
and the fluctuations satisfy S(f{0) = S(p{f3) = . At this stage we mention some technical 



points ||17| , ^8| . First, the part of the action which describes tunneling is equivalent to a one- 
dimensional classical XY model with long range interaction where i3Eq plays the role of 
the system size. Then each update of the phase 5v^(tj) requires a summation over the whole 
discretized lattice. We update Sip{Ti) using the scheme proposed in Ref. [l^ but in a standard 
Metropolis algorithm. Second, in order to extract E^ we need 1 -C f3E^ ~ (3E^ e~°'° . 
Hence large values of do require large ^E^ i.e. a large number of lattice sites. Third, it 
is apparent from Eq. (12) that Z{m) is proportional to the probability for the path to 
visit the m-th sector. Therefore we cannot perform separate simulations for each sector (we 
would trivially obtain Z{m) = 1). On the other hand simulating continuous and discrete 
variables altogether is a technically difficult task. The most economical scheme we found is 
to calculate 

8 



Z{m) JVSipe''^"' Q'Sm-Sm.+i 

sector by sector. Notice that we are interested to the marginal probabihty for a jump to a 
different m sector whereas in the Monte Carlo procedure we update {5ip,m) -^ {6(f',m'). 
This move can be carried out in several ways, according to how we choose the final 6(p' 
configuration. Also the relaxation of the path 6ip can take place by updating sequentially 
or randomly in the lattice. The results we present do not depend on the details described 
above. In each simulation we measure 300 sample points per winding number. Between two 
sampling points the system had time to evolve for 7 more sweeps. We measured the Gibbs 
energy at lower and lower temperatures, i.e. increasing the number of lattice sites for the 
path 6(f up to 1000 at the lowest temperatures. One can obtain the normalized Z{m) from 
the results Eq. (13) by recursion relations and extract the ground state energy Eg{Qx) from 
finite size scaling. Further details are given in [|1^. 

Reliable results are obtained close to Qx ~ , and the effective bandwidth (inverse 
capacitance) has been extracted from the curvature of the band (see Figs. 4 and 5). However 
the region close to the band edges is beyond the capability of the Monte Carlo method. 
Indeed 0{1) terms in the winding number summation Eq.(12) enter with the oscillating 
phase factor exp {27iimQx/e) . Close to the bottom of the band Qx ^ , and the oscillations 
are slow, whereas for Qx ~ e/2 they are fast. At the band edges one expects Z{Qx ~ e/2) ~ 
exp {—(3 {Eq — A)) ^ 1 . Hence one has to extract an extremely small number by adding 
many 0{1) terms with oscillating phase factors. This makes the numerical procedure very 
unstable near Q^ = e/2 . 

We stress that the same argument applies to any analytic calculation which attempts to 
estimate 2(m,) and perform the winding number summation. Examples will be discussed 
below. 



Panyukov and Zaikin |TT| studied the strong coupling regime by a non-standard instanton 
technique in the phase representation Eqs. (2-3). They obtained a renormalized bandwidth 
Eq ~ Ec o? exp (—a) , which agrees with the results of large a scaling, as well as the Monte 



Carlo data down to surprisingly low values of ao • The result is shown in Fig. 5 where the 
curve of Panyukov et al. has been rescaled by a factor ~ 2, i.e. also the accuracy in the 
pre-exponent is remarkable. 

Panyukov and Zaikin further concluded that the ground state energy completely flattens 
in a wide interval around Q^ = ±e/2 . This latter conclusion differs from our picture and 
below we will argue that the possible reason for this discrepancy is the lack of accuracy of 
the instanton calculation which becomes crucial near the degeneracy points. 



The technique presented in [0] is unusual in two senses: i) each instanton possesses 
the standard zero mode related to the location, but also a non-standard one related to their 
width. Thus instantons of all lengths enter Z{Qx) with similar weight, ii) A summation over 
winding numbers m connected with oscillating phase factors has to be performed (cf. Eq. 
(12)). Hence, near the band edges innocent looking approximations can profoundly alter 



the result. We reinvestigated the accuracy and consistency of the method used in Ref. [jTT 



In particular a dilute instanton approximation has been used which in standard fixed size 
instantons calculations is justified when the A/"- instanton configurations playing the main 
role have finite M . Then for fixed size instantons the dilute limit is always reached for 
/3 ^ oo . As the size of the instantons is not limited here, we have to calculate the average 
width {(j)jsf of an instanton in an AT-instanton configuration and the dilute limit is reached 
only ii M{(j)j\f ^ (3 for the relevant values of A/". In the partition function for A/'-instantons 
Z{M) we first separate the zero modes (the locations Tj and the width ctj of each instanton) 
and integrate the remaining fluctuation determinant to get 

<"->^ = ^ Z{W) J n *^ / n "a, a. = ^^^ . (14) 



As in Ref. [^ we constrain the integration to non-overlapping instanton configurations 
{tj + o-j/2 < Tj^i — o-j+i/2). Hence the result (14) is a lower limit to the average size 
and M{(j)j^ is never small compared to (3 . We conclude that the instanton gas is never 
dilute. Thus the contributions from overlapping instanton configurations are comparable to 
the non-overlapping ones considered in [lTl| . 
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We calculated also the bare interaction between an instanton and an anti-instanton 
S{i,j) = — 16a;cr.jcrj/[(crj + (Tj)^ + (tj — Tj)^] which then attract each other. Instantons 
of the same "sign" do not interact. Close pairs of instantons are favoured, and this will 
amplify the deviations from the dilute, non-interacting instanton gas picture. Since close to 
the degeneracy points extreme accuracy is required we consider the result at the edges of the 



bands derived in Ref. |]Tl|], which are based on considerable approximations, not conclusive. 
However close to the bottom of the band results are less sensitive to approximations, which 
explains the excellent agreement between the instanton calculation and our Monte Carlo 
results. 

The exponential suppression of the bandwidth is also found in the straight-line approx- 



imation, introduced in |T^. We make use of the decomposition of v'(t), introduced for the 
Monte Carlo simulation, and expand the action in terms of the fluctuations 5^5 (r) around 
the straight lines 2TTmT/f3 + (fo 



2 ~ 

Sm{S(p) ^ a |m| + T V Iytt + 7^ {\^u - ^ml + \^u + ^m\ - 2|t^m|)] \6Lpy\^ . (15) 

This action gives rise to soft modes for IcUj^l < \ujm\- The fluctuation determinant relative to 
the m = term is 

where a = ^^^- In the a ^ 1 limit, we obtain the partition function 

m 

In the temperature regime Ec<ye~°' ^ T ^ Ec we can confine ourselves to m = 0, ±1. 
This yields (1/(3) \nZ{Qj.) k, const + E^ cos {liiQ^je) where E*q = ^^^ e~" in agreement 
with the previous results. However at lower temperatures the approximation breaks down 
and the numerical summation over m gives a negative partition function. This breakdown is 
understandable since it happens at the crossover temperature where according to the scaling 
analysis the effective a decreases below 1, and quantum fluctuations of Sip{T) become large. 
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V. OTHER RESULTS AND DISCUSSION 

A scaling analysis, very similar to the one we present in section III, was performed in Ref. 
0. The gap renormalization of Eq. (7a) had been found there, but the possibility of wave 
function or a renormalization was not considered. As a result a phase transition between a 
finite gap and a zero gap region at a^ = 1/2 was predicted. The main consequence of the 
additional scaling of a of Eq.(7b) is that the gap remains finite. The transition is smeared, 
leaving only a strong crossover around ao ~ 1/2 (see Fig. 4). 

Two studies addressed directly the weak coupling regime, one performing a poor man's 
scahng analysis JlOl, and one solving a Dyson equation |^. Our results agree with those in 



the leading logarithmic approximation. Difference arises because in our formalism the gap 
appears as an explicit low energy cutoff, yielding a self-consistent equation for A. This will 
be important when Oq is not very small. 

The ground state energy Eq, the voltage at the junction C{V) = dEG{Qx)/dQx and 
the average number (n) of electrons on the island are shown in Figs. 6, 7 and 8. As Ref. 



Ill found a complete flattening of the band around the degeneracy points, it was suggested 
that in the middle of the vertical part of the original sawtooth pattern of (V^) pO[ a new 
S'-shape develops. In the light of the above analysis there is no support for the complete 
flattening. As no reliable treatment of the gap is available in the large a regime, we used 
the scaling formalism in the ao ~ 1 regime (see Fig. 7). The logarithmic corrections of 
Eq.(7) modify the vertical part in a weak manner. The clearest consequence is the strong 
suppression of the amplitude of the sawtooth oscillations already at moderate values of ao- 
Another observable feature is the possible broadening of the excited states, caused by 
the appearance of a finite lifetime F"^. We adopt the method of Ref. to determine F 
close to the band edges, in the ao < 1 regime. Recall that for Q^ < e/2 the ground state 
energy close to the degeneracy point is given by EciQx) ~ -Ec/4 — A{Qx)/2. We can 
determine the energy of the first excited state Ei{Qx) by analytic continuation of Ec{Qx) 
to Qx > e/2 . The argument of the logarithm in A, turns negative when passing Qx = e/2 
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and an imaginary part develops in Ei{Qx)- The inverse lifetime is then determined by 

r = 9™(£,wj)«- *-^^^, . (18) 

' 'l + 2a„ hi(Bc/A„) + (4to„)2 



As r ^ A, the excited levels remain well defined [21|. The primary experimental conse- 
quence of r > is that the I-V curves become smoothed proportional to F/A around the 
onset of the Coulomb-blockade. 

At finite temperatures also thermal fiuctuations have to be considered. If T <^ Ec and 
for ao < 1 the physics near the band edges still involves only the two lowest charge states 
and the Gibbs energy is given by 

G(T, Q,) ^ EoiQ.) + Ao/2 - -^ In [ 2 cosh (^)] (19) 

The renormalized gap is then calculated using Eq.(8) with Uc = T for T > A{Qx) 
and with c^c = A for T < A{Q^) . All the quantities of interest can be calculated and 
in particular the normalized {V) vs. Qx has now a finite slope at Qx = e/2 , given by 
{^PEc)/[l + 2q;o ln(/5i?c)]^ • This result has also been found by other methods |^. Notice 
that a finite tunneling strength is very effective in suppressing the slope when 2 ao \n{i3Ec) ~ 
1 . In the experiments of the Saclay group [Q, 2ao ln(/?-E'c) ~ 10^^ , so quantum fiuctuations 
of the charge do not explain the observed suppression of the slope at Qx = e/2. In this case 
the screening of the charge is probably due to the fact that thermal noise coming from the 
electromagnetic environment can excite tunneling |2^ . 

Finally we reconsider the case of S-S junctions with finite subgap quasiparticle tunneling. 
Here the effect of a Josephson coupling between the electrodes has to be considered as well. 
For large Ej a Kosterlitz-Thouless type transition was found in Ref. [^: for ao > 1/4 
the Josephson coupling Ej scales towards larger values. Here we discuss the small Ej 
limit, where Ej does not renormalize (see below). The flow diagram in the Ej/Eq-cx is 
shown in Fig. 9. plane. In the strong tunneling regime the flow lines of are given by 
Ej/ Ec oc exp(— CKo) which is due to the exponential suppression of Ec discussed in 
section IV. One of these lines is the separatrix between the two phases which ends at a = 1/4 
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for Ej -^ oo. In the regime a < 1 (i.e. a < 27r^) we can study at the charge representation 
Eq. (4), modified by the effect of the Josephson tunnehng @J^. The typical configurations 
differ from those shown in Fig. 3 because 2e bhps and antibhps are present, due to the 
transfer of Cooper pairs. However, no infrared process is connected with the transfer of 
Cooper pairs and hence Ej does not renormalize. 

We thus arrive at the following picture for the the small Ej regime: for large a the 
system flows towards larger Ej/ Ec until it reaches a < 1 (a < l/27r^) where Ec itself does 
not scale appreciably anymore, and the flow lines flatten. In general a will scale towards 
smaller (cf. Eq.(8)) but flnite values, as discussed in section III. Thus the flow lines in Fig. 
9 do not reach the a = axis. The actual flnal value of a depends on the single particle gap 
^q{Qx)- Also the detailed behaviour of the separatrix in the intermediate a and Ej regimes 
may depend on Aq. 

We acknowledge discussions with E. Brezin, R. Fazio, D. Esteve, R. Fresard, H. Grabert, 
F. Guinea, H. Pothier, H. Schoeller and A.D. Zaikin. This work is part of the Sonder- 
forschungsbereich 195, supported by the DFG. GTZ has been supported by NSF grant 
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Figure Captions : 

Fig. 1 : The single electron box. It shows the interplay between the discrete charge —ne 
on the island and the continuous charge Qr^ = CgVx controlled by the voltage source. 

Fig. 2 : The band structure for ao = 0. The energy gap Aq can be tuned between Ec and 
zero by varying the external voltage. 

Fig. 3 : Some diagrams contributing to Z. Full lines represent the allowed trajectories in 
the charge space, wiggly lines are associated with a{Tj — rj). Repeated blip-antiblip pairs 
(a), rainbows (b) and crossings (c) are retained in for the TLS approximation. The diagram 
in (d) involves four charge states. 

Fig. 4 : The effective gap A and tunneling strength a obtained by integrating the scaling 
equations for different initial values {Aq/Eq = 0.01 , 0.05 , 0.1 ), compared with the values 
inferred from Monte Carlo simulations (diamonds). 

Fig. 5 : Monte Carlo results (diamonds) for the effective capacitance renormalization in 
the weak tunneling regime and in the (non-perturbative) intermediate tunneling regime. 



Comparison is made with the results from perturbation theory, from Ref. [|TT| and from Ref. 
10. 

Fig. 6 : The renormalized energy bands close to the edges for various ao in the weak coupling 
limit. 

Fig. 7 : The normalized voltage at the junction C{V)/e = {Q)/e close to the band edges, 
at T = , for various cto in the weak coupling limit. 

Fig. 8 : The expectation value of the number of excess electrons in the island close to the 
band edges, at T = , for various Oq in the weak coupling limit. In the absence of charging 
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effects the "ohmic" linear dependence is found. 

Fig. 9 : Flow diagram in the ^ - a plane for S-S tunnel junctions with finite subgap 
conductance ao- 



